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Abstract. We use a measure of clustering derived from 
the nearest neighbour distribution and the void probabil- 
ity function to distinguish between regular and clustered 
structures. This measure offers a succinct way to incor- 
porate additional information beyond the two-point cor- 
relation function. Application to a supercluster catalogue 
by Einasto ct al. (1997d) reveals no clustering in the dis- 
tribution of superclusters. However, we show that this su- 
percluster catalogue is severely affected by construction 
effects. Taking these biases into account we still find some 
indications for regularity on the largest scales, but the sig- 
nificance is only at the one-cr level. 

Key words: large-scale structure of the Universe - Cos- 
mology: theory - Galaxies: clusters: general 



1. Introduction 



In a recent paper Einasto et al. (1997c) report a peak in 
the 3D-power spectrum of a catalogue of clusters on scales 
of 120/i _1 Mpc. Broadhurst et al. (1990) observed period- 
icity on approximately the same scales in an analysis of 
ID-data from a pencil-beam redshift survey. As is well 
known from the theory of fluids, the regular distribution 
(e.g. of molecules in a hard-core fluid) reveals itself in 
an oscillating two-point correlation function and a peak 



in the structure function respectively (see e.g. Hansen & 



McDonnald 1986). In accordance with this an oscillating 



two-point correlation function £2 (?") or at least a first peak 



was reported on approximately the same scale ( Kopylov et 



al. 1988, Mo et al. 1992, Fetisova et al. 1993, and Einasto 



et al. 1997a|) 



In this paper we analyze the supercluster catalogue of 



Einasto ct al. (1997d) which was constructed from an ear- 



lier version of the cluster catalogue by Andernach & Tago 



(1998) using a friend-of-friends procedure. With methods 



based on the nearest neighbour distribution and the spher- 
ical contact distribution we can show that this superclus- 
ter catalogue is regular with 95% significance. However, 
taking into account the selection and construction effects, 
the high significance vanishes and we only find some indi- 



cation for a regular distribution on large scales, showing 
that this supercluster catalogue is seriously affected by the 
construction method with a friend-of-friends procedure. 

This paper is organized as follows. In Sect. || we dis- 
cuss our methods. Through some examples, we illustrate 
the properties of the J-function and show that it offers a 
concise way to incorporate information about correlations 
of arbitrary order. The analysis of the supercluster dis- 
tribution and of a set of mock supercluster catalogues is 
presented in Sect. |[ We summarize our results in Sect. [l| 

2. Methods 

To analyze the set of points X = {xi},^, 6 K 3 given 
by the redshift coordinates^ of the superclusters we use 
the spherical contact distribution F(r), i.e. the distribu- 
tion function of the distances r between an arbitrary point 
and the nearest object in X. F(r) is equal to the expected 
fraction of volume occupied by points which are not far- 
ther away than r from the objects in X. Therefore, F(r) 
is equal to the volume density of the first Minkowski func- 



tional as introduced into cosmology by Mecke et al. (1994) 



As another tool we use the nearest neighbour distribution 
G(r), that is defined as the distribution function of dis- 
tances r of an object in X to the nearest other object in X. 
For a Poisson distribution the probability to find a point 
only depends on the mean number density /?, leading to 
the well-known result 



Gp (r) = 1 — exp 



_47T , 

H 3 



Mr)- 



(1) 



Recently, van Lieshout & Baddeley (1996) suggested to 
use the ratio 



J(r) = 



1 - G(r) 
1 - F(r) 



(2) 



as a probe for clustering of a point distribution. For a Pois- 
son distribution J(r) = 1 follows directly from Eq. (|I|). 
As shown by van Lieshout & Baddeley (1996), a clustered 



1 Throughout this article we measure length in units of 
/i _1 Mpc, with Ho = WOh km s _1 Mpc" 1 . 
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point distribution implies J(r) < 1, whereas regular struc- 
tures are indicated by J{r) > 1. Typical clustered struc- 



tures are produced by Neyman-Scott processes (Neyman 



& Scott 1958) which have been used to model the distri- 
bution of galaxies. Regular structures arc seen in a peri- 
odic, or a crystalline arrangement of points. In a statisti- 
cal sense, and opposed to clustering, regular ("ordered") 
structures are also seen in liquids. Qualitatively one may 
explain the behaviour of J(r) as follows: 

— In a clustered distribution of points G{r) increases 
faster than for a random distribution of points, since 
the nearest neighbour is typically in the close sur- 
roundings. F(r) increases more slowly than for a ran- 
dom distribution, since an arbitrary point is typically 
inbetween the clusters. These two effects give rise to 
J(r) < 1. 

— On the other hand, in a regular distribution of points, 
G(r) increases more slowly than for a random distribu- 
tion of points, since the nearest neighbour is typically 
at a finite characteristic distance (e.g. in the case of 
a crystal). F(r) increases faster, since the typical dis- 
tance from a random point to a point on a regular 
structure is smaller. These two effects cause J(r) to be 
greater than unity. 

— J(r) = 1 indicates the borderline between clustered 
and regular structures. 

2.1. Gaussian approximation 

For a stationary point distribution 1 — F{r) is equal to 
the probability that no galaxy is inside a sphere B r with 
radius r: 



F(r) = l-p (B r ), 



(3) 



with po{B r ) being the void probability function ( White] 



197S ). Similarly, 1 — G{r) is equal to the probability that 



there is an object at x e X and that there is no other 
object inside the sphere B r {n) centered on x: 



G{r) 



pi(x | B r (s)) 



(4) 



/ oi(x | £> r (x)) is the density that we observe a point at x 
under the condition that B r (x) is empty^j. Therefore we 
obtain ( (Sharp 198l| ): 



J(r) 



Pi(x | g r (x)) 
P 



(5) 



Following Stratonovich (1963) we can express the condi- 
tional density p\ in terms of the n-point correlation func- 



2 Assuming a stationary point distribution we can choose x 
to be the origin. 
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Fig. 1. Above we see the two-point correlation function 
£2(7") of an anticorrelated point distribution (dashed), a 
correlated point distribution (dotted), and of a process 
mimicking the cluster distribution (solid); r in arbirary 
units. Below the corresponding J(r) with p = 1 in the 
Gaussian approximation are depicted. 



tions, the normed cumulants, (see also |White 1979| ): 
Pl (x\B r ( X ))=p-J2- ' ' 



.1! 

n=l 

x / d 3 xi--- d 3 x n £ n+ i(x,xi, . . . ,x„). 

'Br(x) JB r (x) 

A Gaussian approximation, i.e. £„ = for n > 2, yields 



(G) 



J(r) 



1 



2vr f ds s 2 6(s)- (7) 
Jo 



In Fig. [I] we show several two-point correlation 
functions £2(^)1 satisfying the normalization condition 
J °°ds s 2 ^(s) = 0, and the corresponding J(r) in the 
Gaussian approximation. From the third example (solid 
lines in Fig. [I]) we see that a point distribution which is 
correlated (clustered) on small scales may show anticor- 
relation on large scales, such that J(r) < 1 for small and 
J(r) > 1 for large r. Two simple examples of a correlated 
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distribution with £2 (r) > 0, and of an anticorrelated distri- 
bution with £2(7") < clearly show the expected J(r) < 1 
and J(r) > 1 respectively. 

2.2. Beyond the Gaussian approximation 

According to Eqs. (||) and @ J(r) depends on correla- 
tions of arbitrary order. Therefore, a Gaussian approxi- 
mation to J(r) may be misleading. We illustrate this with 
two different point distributions: a Poisson (i.e. random) 
distribution of points, and points given by the model of 
Baddeley & Silverman (1984). Both point distributions ex- 
hibit the same two-point characteristics, but the example 
of Baddeley fc Silverman (1984)| is regular by construc- 
tio n. To generate a realization o f the point distributions 
by Baddeley fc Silverman (1984) we divide the unit square 
into 20 x 20 cells and randomly place 0, 1, or 10 points 
into each cell, with a probability of 1/10, 8/9, and 1/90 
respectively. In Fig. ^ we display these point distributions 
with 388 points in a square. By visual inspecti on the set of 
points given by Baddeley fc Silverman (1984) shows a reg- 
ular structure. Larger voids are only seen in the Poisson 
distribution. We obtain J(r) > 1 for the regular point set, 
clearly distinguishable from the J(r) = 1 for the Poisson 
distributed points (Fig. |J). Both have the same two point 
correlation function ^z{r) = by construction. No scale 
can be deduced from the two-point correlation function 
£2 (f) as seen in Fig. |[ Since the number density and £2 (?") 
are equal in both point sets, the difference in J(r) results 
from high-order correlations only. J(r) and its variance 
diverge near the intrinsic scale r = 1/20 of this specific 
regular point distribution, since F(r) approaches unity. 

If we estimate F(r) from one realization of a point 
process, F(r) becomes unity when r becomes larger than 
the radius of the biggest empty sphere which fits inside 
the sample. Similarly, we get G(r) = 1 for r larger than 
the largest distance between neighbouring objects. There- 
fore the estimate of the J-function from a single point 
set becomes undetermined beyond these radii. Since J(r) 
is a global measure we are still able to detect regular 
structures, as global features of a point pattern, from the 
J(r) > 1 for radii r below the intrinsic scale, here 1/20. 



As discussed in Bedford fc van den Berg (1997) , J(r) 
1 does not necessarily imply Poisson distributed points. 



The morphological measure J(r) was used by Kerscher 



et al. (1997) to investigate large scale fluctuations in the 
galaxy distribution and earlier by [Sharp (1981) 



to test a 



hierarchical model for galaxy clustering. 

Since all real astronomical catalogues are spatially 
limited we have to use edge-corrected estimators as de- 
tailed in Kerscher et al. (1997). The rationale behind 



these estimators is to use only the points whose possible 
nearest neighbours are contained in the sample window. 
With these estimators we neither make any assumptions 
about the exterior of our sample, nor do we use weighting 



;/'■ (a) 



:(b) 



Fig. 2. In plot (a) randomly distributed points and in 
plot (b) points from realizatio n of the regular example by 
Baddeley fc Silverman (1984) are displayed. 



schemes. However, our analysis of the supercluster sample 
is restricted to a radial distance of at most 60-70/i _1 Mpc. 

3. Results 

3.1. The supercluster sample 



We use the supercluster sample compiled by Einasto et 
|al. (1997d)| with 220 superclusters out to z = 0.12. The 
sample was generated with a friend-of-friends algorithm 
using a linking length of 24/i~ 1 Mpc from an earlier version 
of the cluster sample by Andcrnach fc Tago (1998)| and in- 



cludes all superclusters of at least two member clusters. A 



detailed discussion of the sample is given in Einasto et al 
|(1997d) . We limit our analysis to a region within galactic 
latitude |6| > 20° and a maximum radial distance of 330 
/i _1 Mpc. As directly suggested by the sample geometry 
we perform our analysis separately for the northern and 
southern parts (in galactic coordinates). 95 superclusters 
enter into the northern part (mainly Abell sample) and 
116 superclusters into the southern part (mainly ACO 
sample). The selection effects are modeled with an inde- 
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Fig. 3. The plots show the two-point correlation function 
£2(7") on top, and J(r) at the bottom for the random points 
(dark shaded) and for the regular example by Baddcley & 
Silvern ian (1984)| (light shaded). The areas correspond to 
the lcr-error estimated from fifty realizations. The dashed 
lines display the J(r) for the point distributions shown in 
Fig. |^; r is in units of the side length of the square. 



pendent radial and angular selection function (see Einasto 



et al. 1997a). 



3.2. Regular structures? 

In Fig. H the values of J{r) for the supercluster distribu- 
tion are plotted together with the average and lcr-error 
of 99 realizations of a (pure) Poisson distribution with 
the same sample size and geometry. A J(r) larger than 
unity, as expected for a regular distribution of the points 
(see Sect. ||), is clearly seen. J(r) for both parts is above 
one, lying outside the la-range of the Poisson distribu- 
tion on scales from 15/i _1 Mpc to 50/i _1 Mpc. The kink 
in J(r) at r — 45/i~ 1 Mpc in the southern part and at 
55ft.~ 1 Mpc in the northern part indicates the typical scale 
on which the nearest supercluster is situated. This agrees 
with the median distance to the nearest poor supercluster 




20 40 

r [hr 1 Mpc] 



Fig. 4. J(r) for the northern part (dashed), the southern 
part (solid), and a Poisson sample (shaded area) are dis- 
played. The curves are smoothed with a triangular kernel 
with a total width of 3/i _1 Mpc. 



of 45/i _1 Mpc as estimated by Einasto ct al. (1997d). As 
discussed in Sect. || J(r) becomes unreliable on scales be- 
yond 60ft. -1 Mpc (see also Fig. ||). With a (nonpar ametric) 
Monte-Carlo test, as described by Besag fc Digglc (1977)| , 
we conclude with 95% confidence that the superclusters 
given by Einasto et al. (1997d)| are not compatible with 
Poisson distributed points with the same number density. 
We will see that this is not decisive, since up to now we 
did not include selection and construction effects. 

To test the influence of the selection effects and the 
construction process we generate 99 "mock superclus- 
ter catalogues" . We start with Poisson distributed points 
within a sphere of 370/i _1 Mpc incorporating the radial 
and angular selec tion effects of the galax y cluster cat- 
alogue given by Einasto et al. (1997a) ; 



then 



we ap- 



ply a friend-of-friends procedure with linking length of 
24/i _1 Mpc to identify the mock "superclusters". As seen 
in Fig. ^, using a friend-of-friends algorithm, we generate 
an empty sphere with radius of at least 24/i _1 Mpc around 
each supercluster center, introducing an artificial anticor- 
relation, leading to J(r) > 1. Therefore, the regularity 
seen in the northern part of the supercluster sample up 
to scales of 60/i~ 1 Mpc is at least partly an artifact of the 
construction. Still the southern part shows a J(r) above 
the mean J(r) of the mock superclusters, mostly outside 
the lc-range, but a definite statement with a significance 
of 95% (roughly 2a) is no longer possible. 
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Fig. 5. The lower solid line is the result for the north- 
ern part, the upper solid line the result for the southern 
part and the dashed line the result for the whole sample. 
The shaded area marks the la-range estimated from 99 
mock supercluster realizations. The inset plot shows the 
difference A between the mean values of J(r) for the pure 
Poisson process and the mock superclusters samples. All 
curves are smoothed with a triangular kernel with a total 
width of 3/i _1 Mpc. 



4. Discussion and Conclusion 

We have shown that the statistica l properties of the su- 
percluster distribution as given by Einasto ct al. (1997d) 
are seriously affected by the construction with a friend-of- 
fricnds procedure. This is not astonishing since the linking 
length of 24/i _1 Mpc is already one fifth of the claimed reg- 
ularity scale of 120/i _1 Mpc. The distinction between reg- 
ular and clustered point patterns with the J-function is 
unambiguous for a homogeneous and isotropic point distri- 
bution. In such a case the borderline is given by J(r) = 1. 
Our procedure for generating mock supercluster samples, 
where we start with a Poisson sample, include the selec- 
tion effects, and redo the supercluster identification with 
a friend-of-friends algorithm in the same way as for the 
real cluster sample, results in a J(r) ^ 1 even though we 
started from a Poisson distribution. Although it is plau- 
sible, that such mock-supercluster samples describe the 
borderline between clustered and regular structures, there 
is no proof of this assertion. The apparent regularity in the 



northc 'n part of the sample can be explained as a result 



The results for the oscillating two-point correlation 
function of galaxy clusters and, correspondingly the peak 
in the power spectrum were obtained with estimators us- 
ing weighting schemes and boundary corrections, which 
rely heavily on the assumption of homogeneity. Up to 
now there is no reliable way to prove this from the three- 
dimensional distribution of galaxies and clusters. There 
are some hints that the universe reaches homogeneity on 
scales a bove several h undre ds of /t~ 1 Mpc (see the di scus- 
sions by Guzzo 1997 versus riylos Labini et al. 1998). We 



adopted a conservative point of view and used estimators 
which do not make any assumptions about the distribu- 
tion of superclusters outside the sample window. In Sect. || 
we showed that the J-function can be estimated from one 
point set only for scales smaller than the radius of the 
largest void. Therefore, we do not reach the claimed reg- 
ularity scale at 120/i _1 Mpc. Still, the measure J(r) gives 
us information about global properties of the superclus- 
ter distribution, in our case a tendency towards regular 
structures. 

We analyzed the distribution of clusters of the more 



recent redshift compilation by Andernach & Tago (1998) 
with the J-function. We found the expected clumping of 
galaxy clusters, as indicated by J(r) < 1. Qualitatively, 
the J(r) < 1 may be explained with a ^(f) > and the 
Gaussian approximation in Eq. (Q). This clumping out 
to scales of 40/i _1 Mpc is confined mainly to the interior 
of the superclusters. Isolated "field" clusters were not in- 
cluded in the supercluster sample but may contribute to 
the correlation seen up to scales of 50/i -1 Mpcin the clus- 
ter samples. One hierarchical level higher, the supercluster 
centers themselves show a tendency towards regular struc- 
tures. Again this can be explained qualitatively with the 
Gaussian approximation (see Fig. [l]) . A theoretical exam- 
ple illustrating such a hierarchical property is given by 
Neyman-Scott processes (Neyman & Scott 1958): In such 
a process the overall distribution of points shows correla- 
tion (i.e. £(r) > for small r), but the cluster centers of 
these points are distributed randomly by construction. 

Unlike the two-point correlation function the J— 
function incorporates information stemming from high or- 
der correlations. Our example in Fig. || illustrates, that 
a regular structure detected unambigously with the J— 
function may not be visible in an analysis with the two- 
point correlation function £2(1") alone. 

Another problem is the fluctuations between the north- 
ern and southern parts of the sample. This may be at- 
tributed to the different selection effects entering the Abcll 
and ACO parts of the sample, probably due to the differ- 
ent sensitivity of the photo plates used. However, in the 
case of the IRAS 1.2 Jy galaxy catalogue such fluctuations 



of these construction effects, whereas the southern part 
still shows a trend towards regular structures outside the 
lcr-range. 



were shown to be real on scales of 200/i 1 Mpc ( Kerscher et 
al. 1997] ). Also > |Zucca et al. (1997)| find from the ESP sur 



vey, that at least in the southern hemisphere the local den- 
sity is below the mean sample density out to 140/i _1 Mpc. 
If we assume that the fluctuations decrease on scales above 
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200/i~ 1 Mpc, the finding of regular structures on such large 
scales is a great challenge to the standard scenarios of 
structure formation by gravitational instability, starting 
from Gaussian initial density fluctuations. Implications of 
these regular structures for the standard scenarii of struc- 
ture formation are discussed in Einasto ct al. (1997b)| and 
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